Multiscale simulations of porous media flows in flow-based coordinate system

In this paper, we propose a multiscale technique for the simulation of porous media flows in a flow-based coordinate system. A flow-based coordinate system allows us to simplify the scale interaction and derive the upscaled equations for purely hyperbolic transport equations. We discuss the applications of the method to two-phase flows in heterogeneous porous media. For two-phase flow simulations, the use of a flow-based coordinate system requires limited global information, such as the solution of single-phase flow. Numerical results show that one can achieve accurate upscaling results using a flow-based coordinate system

Robust Execution of Contact-Rich Motion Plans by Hybrid Force-Velocity Control

In hybrid force-velocity control, the robot can use velocity control in some directions to follow a trajectory, while performing force control in other directions to maintain contacts with the environment regardless of positional errors. We call this way of executing a trajectory hybrid servoing. We propose an algorithm to compute hybrid force-velocity control actions for hybrid servoing. We quantify the robustness of a control action and make trade-offs between different requirements by formulating the control synthesis as optimization problems. Our method can efficiently compute the dimensions, directions and magnitudes of force and velocity controls. We demonstrated by experiments the effectiveness of our method in several contact-rich manipulation tasks. Link to the video: https://youtu.be/KtSNmvwOenM.Comment: Proceedings of IEEE International Conference on Robotics and Automation (ICRA2019

Algebraic Reduction of Feynman Diagrams to Scalar Integrals: a Mathematica implementation of LERG-I

A Mathematica implementation of the program LERG-I is presented that performs the reduction of tensor integrals, encountered in one-loop Feynman diagram calculations, to scalar integrals. The program was originally coded in REDUCE and in that incarnation was applied to a number of problems of physical interest.Comment: 16 page

A modified particle method for semilinear hyperbolic systems with oscillatory solutions

We introduce a modified particle method for semi-linear hyperbolic systems with highly oscillatory solutions. The main feature of this modified particle method is that we do not require different families of characteristics to meet at one point. In the modified particle method, we update the ith component of the solution along its own characteristics, and interpolate the other components of the solution from their own characteristic points to the ith characteristic point. We prove the convergence of the modified particle method essentially independent of the small scale for the variable coefficient Carleman model. The same result also applies to the non-resonant Broadwell model. Numerical evidence suggests that the modified particle method also converges essentially independent of the small scale for the original Broadwell model if a cubic spline interpolation is used

Preconditioning Markov Chain Monte Carlo Simulations Using Coarse-Scale Models

We study the preconditioning of Markov chain Monte Carlo (MCMC) methods using coarse-scale models with applications to subsurface characterization. The purpose of preconditioning is to reduce the fine-scale computational cost and increase the acceptance rate in the MCMC sampling. This goal is achieved by generating Markov chains based on two-stage computations. In the first stage, a new proposal is first tested by the coarse-scale model based on multiscale finite volume methods. The full fine-scale computation will be conducted only if the proposal passes the coarse-scale screening. For more efficient simulations, an approximation of the full fine-scale computation using precomputed multiscale basis functions can also be used. Comparing with the regular MCMC method, the preconditioned MCMC method generates a modified Markov chain by incorporating the coarse-scale information of the problem. The conditions under which the modified Markov chain will converge to the correct posterior distribution are stated in the paper. The validity of these assumptions for our application and the conditions which would guarantee a high acceptance rate are also discussed. We would like to note that coarse-scale models used in the simulations need to be inexpensive but not necessarily very accurate, as our analysis and numerical simulations demonstrate. We present numerical examples for sampling permeability fields using two-point geostatistics. The Karhunen--Loève expansion is used to represent the realizations of the permeability field conditioned to the dynamic data, such as production data, as well as some static data. Our numerical examples show that the acceptance rate can be increased by more than 10 times if MCMC simulations are preconditioned using coarse-scale models

A 3D Numerical Method for Studying Vortex Formation Behind a Moving Plate

In this paper, we introduce a three-dimensional numerical method for computing the wake behind a flat plate advancing perpendicular to the flow. Our numerical method is inspired by the panel method of J. Katz and A. Plotkin [J. Katz and A. Plotkin, Low-speed Aerodynamics, 2001] and the 2D vortex blob method of Krasny [R. Krasny, Lectures in Appl. Math., 28 (1991), pp. 385--402]. The accuracy of the method will be demonstrated by comparing the 3D computation at the center section of a very high aspect ratio plate with the corresponding two-dimensional computation. Furthermore, we compare the numerical results obtained by our 3D numerical method with the corresponding experimental results obtained recently by Ringuette [M. J. Ringuette, Ph.D. Thesis, 2004] in the towing tank. Our numerical results are shown to be in excellent agreement with the experimental results up to the so-called formation time

Flux rope proxies and fan-spine structures in active region NOAA 11897

Employing the high-resolution observations from the Solar Dynamics Observatory (SDO) and the Interface Region Imaging Spectrograph (IRIS), we statistically investigate flux rope proxies in NOAA AR 11897 from 14-Nov-2013 to 19-Nov-2013 and display two fan-spine structures in this AR. For the first time, we detect flux rope proxies of NOAA 11897 for total 30 times in 4 different locations. These flux rope proxies were either tracked in both lower and higher temperature wavelengths or only detected in hot channels. Specially, none of these flux rope proxies was observed to erupt, but just faded away gradually. In addition to these flux rope proxies, we firstly detect a secondary fan-spine structure. It was covered by dome-shaped magnetic fields which belong to a larger fan-spine topology. These new observations imply that considerable amounts of flux ropes can exist in an AR and the complexity of AR magnetic configuration is far beyond our imagination.Comment: 8 pages, 8 figures, Accepted for publication in A&

Second-Order Convergence of a Projection Scheme for the Incompressible Navier–Stokes Equations with Boundaries

A rigorous convergence result is given for a projection scheme for the Navies–Stokes equations in the presence of boundaries. The numerical scheme is based on a finite-difference approximation, and the pressure is chosen so that the computed velocity satisfies a discrete divergence-free condition. This choice for the pressure and the particular way that the discrete divergence is calculated near the boundary permit the error in the pressure to be controlled and the second-order convergence in the velocity and the pressure to the exact solution to be shown. Some simplifications in the calculation of the pressure in the case without boundaries are also discussed

Multiscale Finite Element Methods for Nonlinear Problems and their Applications

In this paper we propose a generalization of multiscale finite element methods (Ms-FEM) to nonlinear problems. We study the convergence of the proposed method for nonlinear elliptic equations and propose an oversampling technique. Numerical examples demonstrate that the over-sampling technique greatly reduces the error. The application of MsFEM to porous media flows is considered. Finally, we describe further generalizations of MsFEM to nonlinear time-dependent equations and discuss the convergence of the method for various kinds of heterogeneities